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ABSTRACT 

Cosmic magnification has been detected through cross correlation between foreground 
and background populations (galaxies or quasars). It has been shown that weighing 
each background o bject by its a — 1 can significantly improve the co smic magnifica- 
tion measurement (jMenard fc Bartelmamill2002t IScranton et al.ll2005fh Here, a is the 
logarithmic slope of the luminosity function of background populations. However, we 
find that this weighting function is optimal only for sparse background populations in 
which intrinsic clustering is negligible with respect to shot noise. We derive the optimal 
weighting function for general case including scale independent and scale dependent 
weights. The optimal weighting function improves the S/N (signal to noise ratio) by 
~ 20% for a BigBOSS-like survey and the improvement can reach a factor of ~ 2 for 
surveys with much denser background populations. 
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1 INTRODUCTION 

Weak gravitational lensing directly probes t he matter distri- 
bution of the universe fe.g. |Refregierll2003l 1) and is emerg- 
ing as one of t he most powerful pr obes of dark mat- 
ter, dar k energy jA lbrccht ct al. 2006) and the nature of 
gravity (|jain fc Zhand 2008). By far the most sophisti- 
cated way to measure weak lensing is through cosmic 
shear, which i s the lensing i nduced coherent distortion in 
galaxy shape (|Fu et al.ll2008l and references therein). Coor- 
dinated projects on precision weak lensing reconstruction 
through galaxy shapes have been car ried out extensively 
(STEP, iHevmans etafl | 2006l; STEP2, iMassev et all 120071; 



Bridle et al 



2009; GREAT10, iKitching et ail 



GRE AT08, 
l2010h . 

Alternatively, one can reconstruct weak lensing through 
cosmic magnification, namely the lensing i nduced coher - 



Here, k is the lensing convergence and a(m,z) = 
2.5 dlog n(m, z) /dm — 1 is a function of the galaxy apparent 
magnitude m and redshift z. The number count of galaxy 
brighter than m is N(m) = J m n(m)dm. Throughout the 
paper we use the superscript "L" to denote the lensed quan- 
tity. 

Since cosmic magnification does not involve galaxy 
shape, weak lensing reconstruction through cosmic magni- 
fication automatically avoids all problems associated with 
galaxy shape. The key step in such reconstruction is to elim- 
inate <5 g , which is often orders of magnitude larger than the 
lensing signal n. Usually this is done by cross correlating 
foreground population (galaxies) and background popula- 
tion (galaxies or qua s ars) with no overlapping in redshift 
dScranton et all 120051 ; iHildebrandt et all 120091 ; IWang et all 
|20H| V The resulting cross correlation is 



ent distortion in galaxy number density ( e .g. Gunnl l|l967t ) ; 
iMenard fc Bartelmannl (120021 1; I Jain et all (|2003h and refer- 
ences therein). Neglecting high order corrections, the lensed 
galaxy (quasar) number overdensity <5g is related to the in- 
trinsic overdensity S s by 

S\ = S s + 2(q - (1) 



2(a b -l)(<5 g ,f(fl)«b(fl')> 



(2) 



Throughout the paper we use the subscript "b" for quanti- 
ties of background population and "f ' for that of foreground 
galaxies. The above relation neglects a term proportional to 
(Kf/tbV, which is typically much smaller than the (<5 g ,fKb) 
termlj 
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This term can be non-negligible or even dominant for suffi- 
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It is important to weigh the cross correlation measure- 
ment appropriately to improve the S/N (signal to noise ra- 
tio). S ince the signal scales as a — 1, iMenard &: Bartelmannl 
(2002) first suggested to maximize the S/N by weighing 
each galaxy with its own a — 1. This weighting significantly 
improves the measurement and a robust 8<r detection of 
the cosmic magnifica tion was achieved for the first time 
jScranton et al.ll2005l ). 

Nevertheless, we find that, the a — I weighting is optimal 
only in the limit where the background galaxy (quasar) in- 
trinsic clustering is negligible with respect to the shot noise 
in background distribution. The statistical errors (noises) 
are contributed by both shot noise and intrinsic clustering 
of foreground and background galaxies. In this letter, we 
derive the exact expression of the weighting function opti- 
mal for the cosmic magnification measurement through cross 
correlation. The new weighting can further improve the S/N 
by ~ 20% for galaxy-quasar cross correlation measurement 
in a BigBOSS-like survey. We can also employ high redshift 
galaxies instead of quasars as background sources which can 
have much larger number density and smaller bias. Smaller 
shot noise results into better performance for the derived 
optimal weighting. The improvement over the a — 1 weight- 
ing can reach a factor of ~ 2 for surveys with background 
population density of ~ 2/arcmin 2 . 

Throughout the paper, we adopt the fiducial cosmology 
as a flat ACDM universe with Qa = 0.728, fi m = 1 — Oa, 
erg = 0.807, h — 0.702 and initial power index n = 0.961, 
which are consistent with WM AP seven years best-fit pa- 
rameters |Komatsu et al.ll201 ll ). 



2 THE OPTIMAL WEIGHTING FUNCTION 

We are seeking for an optimal weighting function linearly op- 
erating on the background galaxy (quasar) number overden- 
sity in flux (magnitude) space. Let's denote the background 
galaxy number overdensity of the i-th magnitude bin as 5^' h 
and the corresponding weighting function as Wi. 

• The simplest weighting function is scale independent, 
so the weighted background galaxy overdensity is 



00 

b 



(3) 



• We can further add scale dependence in Wi to increase 
the S/N. The new weighting function convolves the density 
field. For brevity, we express it in Fourier space as Wi(£). The 
Fourier transformation of the weighted background overden- 
sity is 



(4) 



Here, <5 gi b is the Fourier component of the overdensity <5 gl b- 
The weighting function Wi(£) is real and only depends on 
the amplitude of the wavevector £ = \£\. It guarantees the 
weighted overdensity to be real. 



ciently high redshift foreground galaxy samples jZhang fc Pen! 
l200d) . 



The S/N of the background-foreground galaxy cross- 
correlation depends on the weighting function, so we use 
the subscript "W" to denote the S/N after weighting. The 
overall S/N can be conveniently derived in the Fourier space, 



(-)' 



')w 



CM-g>2 



(2£+l)Al/ 3ky <<^- 



- + (<C gjb >w + (C Sib ) w )(C gif +C aif ) 
(2l+l)Al/ Bky 



l+(C g ,f+C M ) 



<b Bib H-)^c mb + (w^)C sJ 

i(W(c b -l))^Ci 



(5) 



Here, C 



CM-g 



is the cosmic magnification-galaxy cross cor- 
A ~ s is the associated sta- 



relation power spectrum and AC 
tistical error. {■ ■ -)w denotes the weighted average of the 
corresponding quantity. We then have 



CM-g 



) w = <2(a b - l)W)G R 



(6) 



Here, C Kbg is the cross correlation power spectrum between 
background lensing convergence and foreground galaxy over- 
density, (uv) is the averaged product of uv, 



(uv) 



J2 t u(mi)v(mi)N h}i 



(7) 



Here, Nh,i is the number of background galaxies (quasars) in 
the given magnitude bin rrii — Amj/2 < m < rm + Am;/2. 

The S/N scales with /^J and / s k y is the fractional 
sky coverage. C s is the shot noise power spectrum and the 
weighted one is (C E ,b)w = {W 2 }C B ,b- C g is the galaxy power 
spectrum. We adopt a simple bias model for the foreground 
and background galaxies. We then have C g ,j = 6g jC m where 
b gt i is the bias of the i-th magnitude bin and C m is the cor- 
responding matter angular power spectrum. The weighted 
background galaxy power spectrum is 

{C g , h )w = (&g,bW^) 2 C m , b . (8) 



2.1 The scale independent optimal weighting 
function 

The optimal weighing function W can be obtained by vary- 
ing the S/N (Eq. [SJ) with respect to W and maximizing it. 
The derivation is lengthy but straightforward, so we leave 
details in the appendix and only present the final result here. 
The optimal weighting function is of the forir0 



W = (a b - 1) + eb g ,b 



(9) 



where the scale independent constant e is determined by 
Eq. IA2I It is a fixed number for the given redshift bin of 
the given survey. In the limit that shot noise of background 
galaxies overwhelms their intrinsic clustering (C s ,b ~S> C gi b), 
e — >■ 0. In this case, the w e ightin g function a — 1 proposed 
bv lMenard fc Bartelmannl (|2002l ) becomes optimal. 



2 The derived scale independent weighting function implicitly 
assumes no scale dependence in the galaxy bias. In reality, the 
galaxy bias is scale dependent and the application of Eq. [9] is 
limited. The exact optimal weighting function applicable to scale 
dependent bias is given by Eq. 1101 
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2.2 The scale dependent optimal weighting 
function 

The weighting function W (Eq. [9| is optimal under the con- 
dition that W is scale independent. If we relax this require- 
ment and allow for scale dependence in W, we are able to 
maximize the S/N of the cross power spectrum measurement 
at each I bin. Clearly, this further improves the overall S/N. 

In this case, W of different I bins are independent to 
each other. This significantly simplifies the derivation and 
we are now able to obtain an analytical expression for W , 



W(l) = (a h - 1) + 



<(a b - l)6 g , b )C In , b (f)/C S; 



l + (^b>Cm,bW/C,,b 



b s , h {l). (10) 



This form is similar to Eq. (|9]), except that it is now scale de- 
pendent. Here again, in the limit C 3 , b SS> C m ,b ~ C'g.b, W 



a — 1 and we recover the result of iMenard &: Bartelmannl 
(2002). This is indeed the case for SDSS background quasar 
sample. 



2.3 The applicability 

Are the derived weighting functions (Eq. [9] & I10|) directly 
applicable to real surveys? From Eq.O&HSl it seems that we 
need to figure out 6 gjb first. Since & gb = C gib /C mi b and C m , b 
is not directly given by the observation, cosmological priors 
or external measurements (e.g. weak lensing) are required 
to obtain the absolute value of 6 g>b . Hence it seems that 
the applicability of Eq. [9] & [10] is limited by cosmological 
uncertainty. 

However, this is not the case. Eq. llOl shows that, it is the 
combination & gb C m ,b determines W. Since C gj b = fr gib Cm, b 
and «{, are directly observable, Eq. [TO] is determined com- 
pletely by observations. Closer look shows that Eq. [5Jis also 
determined by the combination b 2 h C m ,b, so the correspond- 
ing weighting is determined completely by observations, too. 
Hence the derived optimal weighting functions are indeed 
directly applicable to real surveys. 

For ongoing and planned surveys such as CFHTLS, 
COSMOS, DES, BigBOSS, LSST, SKA, etc., the number 
density of background populations (galaxies) can be high 
and the intrinsic clustering can be non-negligible or even 
dominant comparing to shot noise. In next section we will 
quantify the improvement of the optimal weighting functions 
for a BigBOSS-like survey and briefly discuss implications 
to surveys with even denser background populations. 



3 THE IMPROVEMENT 



BigBOS^fl is a planned spectroscopic redshift survey of 
24000 deg 2 (BigBOSS-North plus South). Cosmic magnifi- 
cation can be measured by BigBOSS through LRG (lumi- 
nous red galaxy)-quasar and ELG (emission line galaxy)- 
quasar cross correlations. In principle, it can also be mea- 
sured through LRG-ELG cross correlation. But the inter- 
pretation of the measured cross correlation signal would 
be complicated by the intricate selection function of ELGs 

3 http://bigboss.lbl.gov/ 



|Zhu et al.ll2009h . In the current paper, we only consider the 
LRG-quasar and ELG-quasar cross correlations. 

There are some uncertainties in the BigBOSS galaxy 
(quasar) redshift evolution, flux distribution and intrinsic 
clustering. To proceed, we will take a number of simplifica- 
tions. So the absolute S/N of cross correlation measurement 
that we calculate is by no means accurate. But our calcula- 
tion should be sufficiently robust to demonstrate the relative 
improvement of the exact optimal weighting function over 
the previous one. 

The LRG and ELG luminosity fu nctions are calculate d 
based on the BigBOSS white paper (|Schlegel et all 120091 ). 
The comoving number density of LRG and ELG is 3.4 x 
W~ 4 (h/Mpc) 3 , then we have 1.1 x 10 7 LRGs in the redshift 
range of z = 0.2 — 1 and 3.3 x 10 7 ELGs in the redshift range 
of z — 0.7 — 1.95. Clustering of LRGs evolves slowly, so we 
adop t LRG bias as & g ,t(z) = 1.7 / D{z) (|Padmanabhan et al.l 
l200rj ). Here D(z) is the linear density growth factor and is 
normalized such that D(z = 0) = 1. Existing knowledge on 
clustering of ELGs i s rath er limited. So we simply follow 
iPadmanabhan et al.l (|2006t ) and approximate the ELG bias 
as 6 g ,f = 0.8/D(z). 

For the luminosity function of background quasars, we 
adopt the LDDE (Luminosity depend ent density ev o lution ) 
model with best fit parameters from ICroom et alj (j2009T) - 
The magnitude limit is g — 23, then we have 2.1 x 10 6 
quasars in the redshift range of z — 2 — 3.5. We choose a 
redshift gap (zb,mm — Zf.max = 0.05) such that the intrinsic 
cross correlation between foreground and background pop- 
ulations can be safely neglected. We adopt a bias model for 
quasar clustering, with &q(^) = 0.53 + 0.289(1 + z) 2 from 
the analysis of 3 x 10 s quasars l|Mvers et alj|2007h . 

The S /N depends on many issues and can vary from 90- 
140 (Table [TJ. The S/Ns of LRG-quasar and ELG-quasar 
correlations are comparable because of a consequence of 
several competing factors including the lensing efficiency, 
galaxy surface density and clustering. Nevertheless, a robust 
conclusion is that BigBOSS can measure cosmic magnifica- 
tion through galaxy-quasar cross correlation measurement 
with high precision. Given such high S/N and accurate red- 
shift available in BigBOSS, it is feasible to directly measure 
the angular di ameter distance from s uch measurement b y 
the method of Ijain fc Taylor! (|2003h : IZhang et all (|2005h : 
iBernsteinl (|20od ). 

Unambiguous improvement in the cross correlation 
measurement by using our optimal weighting (Eq. [9] & I10[) 
is confirmed, as shown in Table [JJ The S/N is improved by 
~ 15% by using the scale independent optimal weighting 
(Eq. [9} and by ~ 20% by using the scale dependent one 
(Eq.HD). 

We further investigate the dependence of the above im- 
provement on the flux dependence of quasar bias. We adopt 
a toy model with b Q (z,F) = b Q (z)(F/F*) fl . Here, F* is 
the flux corresponding to that of the M* in the quasar lu- 
minosity function. /3 is an adjustable parameter and we will 
try P = 0, 0.1, 0.2, 0.3, then the corresponding parameters in 
scale independent weighting are e = 0.049, 0.048, 0.047 and 
0.045. Larger value of ft (^ 0.4) leads to too large quasar bias 
(^ 10) and hence will not be investigated here. Table[T]shows 
consistent improvement by our optimal weighting functions. 
Hence, despite uncertainties in quasar modeling, we expect 
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Table 1. Improving the cosmic magnification measurement by the optimal weighting function. The target survey is BigBOSS. The terms 
on the left side of the arrows are the estimated S/N using the weighting a — I. The ones on the right side are the S/N using the optimal 
weighting function, where the ones on the left side of "/" are what expected using the scale independent weighting (Eq. [9J, and the ones 
on the right are what expected using the scale dependent weighting (Eq. llOll , The improvement depends on the bias dependence on galaxy 
luminosity. To demonstrate such dependence, we adopt a toy model feg oc F@ and investigate different values of the parameter /3. In 
general, the optimal weighting function improves the S/N by 10%-20% for BigBOSS, whose background quasar density is ~ 0.02/arcmin 2 . 
The improvement can reach a factor of ~ 2 for surveys with background (galaxy) populations reaching surface density of ~ 2/arcmin 2 . 



Flux dependence of quasar bias model 


/3 = 


0.1 


0.2 


0.3 


Detection significance of LRG X quasar 
Detection significance of ELG X quasar 


111.2^129.0/136.8 
94.3^106.4/110.3 


109.8-5-126.3/133.9 
93.2->104.5/108.6 


108.3-^123.7/131.1 
92.0-5-102.2/106.7 


106.5^119.3/126.7 
90.5-5-99.1/104.1 



our optimal weighting function to be useful to improve cos- 
mic magnification measurement in the BigBOSS survey. 

Nevertheless the improvement is only moderate. The 
major reason is that, even for BigBOSS, the quasar sample is 
still sparse, with a surface number density ~ 0.02/arcmin 2 . 
Hence shot noise dominates over the intrinsic clustering. In- 
deed, we find that typically C m ,Q/C s ,Q < 0.1. For imaging 
surveys like CFHTLS, COSMOS, DES and LSST, we can 
also use high redshift galaxies as background galaxies to 
correlate with low redshift foreground galaxies. For these 
surveys, high redshift galaxy population (e.g. with z > 1-2) 
can reach surface number density ~ 0.2-2/arcmin 2 or even 
higher. So the shot noise in these surveys can be suppressed 
by a factor of 10-100 or more. The overall improvement of 
our optimal weighting would be larger. 

To demonstrate these further improvements, we hypo- 
thetically increase the surface density of BigBOSS quasars 
by a factor of 10 and 100 respectively, but keep /3 = and all 
other parameters unchanged. Shot noise will be decreased by 
a factor of 10 and 100 correspondingly. The scale indepen- 
dent weighting parameter e can reach 0.12 and 0.22 respec- 
tively. For the first case, the S/N is improved by ~ 38% for 
the scale independent optimal weighting and by ~ 51% for 
the scale dependent one. For the second case, the improve- 
ment is ~ 72% for the scale independent optimal weighting 
and is ~ 94% for the scale dependent one. It is now clear 
that for measuring cosmic magnification through cross corre- 
lation between foreground and background galaxies of many 
existing and planned surveys, one should adopt the optimal 
weighting function derived in this letter. 



4 SUMMARY 

We have derived the optimal weighting functions for cos- 
mic magnification measurement through cross correlation 
between foreground and background populations, for scale 
independent and scale dependent weights respectively. Our 
weighting functions outperform the commonly used weight- 
ing function a — 1 by ~ 20% for a BigBOSS-like survey and 
by larger factors for surveys with denser background popu- 
lations. Hence we recommend to use our optimal weighting 
function for cosmic magnification measurement in BigBOSS, 
CFHTLS, COSMOS, DES, Euclid, LSST, SKA, WFIRST, 
etc. 
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APPENDIX A: DERIVING THE OPTIMAL 
WEIGHTING FUNCTION 

In Section 2, we give the optimal weighting function without 
derivation. Here, we present a brief derivation for the scale 
independent weighting function. The derivation of the scale 
dependent weighting function is similar, but simpler. 
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Maximizing the S/N requires the variation of (S/N) 2 
with respect to W to be zero. From this condition, we obtain 

117 J2 f (2e+l)At v ^(b gtb W) 2 + (W 2 ))/(l+F) 2 

W = — j— '- (a b -l) 

(W(a b -1)> ^(2<+l)A^/ (1+F)2 

<W(a b -l))<b 8-b lV) ^2 e (2t+l)Ae V i,/ (1+F) 2 

7 I 7 ^gib ■ I Al ) 

<VK(a b -l)> ^2 t (2t+l)AiT,/ (1+F) 2 

Here, for brevity, we have denoted ij = C B , b (Cg,f + 
C s . { )/Cl h9 , v = C m , b /C s .b and F (w) = V [u(b e . b W) 2 + 
{W 2 )]/[4{W(a b - I)) 2 ]. 

Noticing that the coefficients of a b — 1 and b gjb only 
involve the weighted average of W and taking the advantage 
that the optimal W remains optimal after a constant (flux 
independent) scaling, we are able to seek for the solution 
of the form W = (a b — 1) + £& g , b , with e a constant to be 
determined. Plugging it into the above equation, we obtain 
the equation of s, 

((a b - l)6g, b > ^ 



(2e+l)Ae v / (2i+l)At V v , h 2 v 

2-<i (i+Fp I L~ii (T+fJ 2 r \ g, b ' 

Noticing that F depends on e. The above equation can be 
solved numerically to obtain the solution of e. 

In the case of scale dependent weighting, each W (£) are 
independent to each other. Through similar derivation, one 
can show that 



w{i) = K - 1) + 



v({a h - l)6 g ,b> 



(A3) 



